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Abstract 

We fitted the vrvr mass distribution in the range 0.5 < M^tt < 1-1 GeV measured in 
hard exclusive positron-proton reactions at HERA by the form dictated by QCD at 
leading twist level. Extracted parameters are related to valence quark distribution 
in the pion, and to the pion and p meson distribution amplitudes. We obtain, 
for the first time, a measurement of the second Gegenbauer coefficient of the p 
meson distribution amplitude: 02"^ = —0.10 it 0.20 for a photon virtuality of (Q^) = 
21.2 GeV^. 

1. Introduction 

Owing to QCD factorisation theorem for hard exclusive reactions []T| the dependence of 
the amplitude of the reaction Q 

72(g) +T(p)^7r+7r-+T'(p') (1) 

with longitudinally polarised virtual photon on the di-pion mass M^^r factors out in an 
universal (independent of the target) factor. At leading order in as(<5^) this factor has 
the form |0]: 

A{M^^)^ / -<l>^(z,C,M^^;g2) . (2) 
Jo z 

Here ^\z, (, M^t^; Q^) is the two-pion light cone distribution amplitude (27rDA) which 
depends on z-longitudinal momentum carried by the quark, ( characterising the distri- 
bution of longitudinal momentum between the two pionsQ, and the invariant mass of 
produced pions M^.^^, the superscript / standing for isospin of produced pions (J = 0, 1). 

iWith Q2 _ 2{p ■ q) > A^cD and M^, Mf, {p' - pf, M^, < . 
^For detailed definition of kincmatical variables z and C see refs. ||. 



The dependence on the virtuahty of the incident photon is governed by the ERBL 
evolution equation 

For the process (|I]) at small Xsj = Q'^/'^ip ■ q) the production of two pions in the 
isoscalar channel is strongly suppressed relative to the isovector channel [0], because the 
former is mediated by C-parity odd exchange. At asymptotically large QCD predicts 
the following simple form for the isovector 27rDA 0, 0: 

$^J„(z,C,M^^) = 6zil - zm - 1)F^(M^.) , (3) 

where Ft^^Mt^t^) is the pion electro-magnetic (e.m.) form factor in time-like region, mea- 
sured with high precision in low energy experiments pO| , |lT|. From eqs. we conclude 
that at asymptotically large QCD predicts unambiguously the shape of the di-pion 
mass distribution. Asymptotically the dependence of the amplitude on M^^r has the form: 

Aasym OC (3{M^^) F^{M^^) COS ^ , (4) 



where /?(M^^) = ^1 — |^ is the velocity of pions in their centre of mass system (cms) 

and 6 is the angle between the directions of the positive pion and the momentum of 
produced tt+tt^ system in the tt+tt" cms. This angle is related to ( in the following way: 



The corresponding di-pion mass distribution has asymptotically the form 

(X M^^ PiM^^)'\F^{M^^)\^ . (6) 



The asymptotic shape for any di-meson production (mesons Mi, M2) was derived in 
0, H, 1^, where it was related to the cross section of e'^e^ Mi, M2 at low y/s. 

At non-asymptotic values, the 27rDA deviates from its asymptotic form (|^). This 
deviation can be described by a few parameters which can be related to quark distributions 
(skewed and usual) in the pion and to distribution amplitudes of mesons (tt, p, etc... ), 
for details see 0. 

2. Deviations from the asymptotic form 

The first non-trivial deviation from the asymptotic form of the 27rDA occurs in P and 
F waves 0. Generically their effect on M^^r dependence of the hard amplitude can be 
written as: 



A{M^^) (X (3{M^^) e'''\F^{M^^)\[l + Di{M^^)j PiicosO) 

+ (3{M^^f e''W2{M^^) Psicos 6) , (7) 

where Pi(cos^^) are Legendre polynomials and 6i{M.„t,) and 63{M.„.„) are the P-wave and 
the F-wave utt scattering phase shifts, which are well known from low-energy experiments. 
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The functions Z^i 2(M7r^) describe the deviation of the amphtude's Mt^t, dependence from 
the asymptotic form. These functions are real and can be parametrised as: 

TTTT 

The dependence of v4i 2((5^) on is governed by the QCD evolution and in leading order 
is given by: 

(OL fnaU 50/(99-6n;) 
^,2(Q^)=^,2(/.0)(^^) . (8) 

With increasing Q^, the parameters Ax ^iQ'^^ go logarithmically to zero and one reproduces 
the asymptotic formula (^. The parameters Ai 2(Q^) are directly related to partonic 
structure of vr and p mesons, see section 4 and ref 0. 

The parameter 63 is Q^-independent but is M^^ dependent. The latter dependence is 
fixed by txtx scattering phase shifts (see for derivation 0): 

h{M^^) = h + Re^ / ds ^. . (9) 

In above equations hi are subtraction constants of corresponding dispersion relations for 
functions Di^2(^7r7r), see details in p|. 

Using expression (0) we can derive the form of the Mtttt distribution: 



dN{M^^^ 
dM^^ 



^ M^^ (3{M^^y D^iM^^Y + (10) 



higher waves / > 5 

where the higher partial waves can be safely neglected. 
3. Angular distributions of produced pions 

Another way to obtain sensitivity to the non-asymptotic parameters Ai^2{Q'^) is to study 
the two pion angular distributions. From the expression of the amplitude (|^) one can 
derive the M^j^ dependence of intensity densities {¥{^{9,4))) defined as: 

d /-^T 



dM„ 



rf(cos0) / d^Yr{e,^)\A{M^^,e,^)\\ [11] 

-1 Jo 



where Yi^{6, if) are spherical harmonics, ip is the azimuthal angle between the pion decay 
plane and plane formed by momentum of the virtual photon and total momentum of 
produced tt+tt^ system. 
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For the first nontrivial intensity density {Y2) we have: 



d 



dM^ 
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F,(M^,)| l + A(A/4,) D2(M,.) cos 



The combination (Fq*^ — \/5/2Y2) is especially sensitive to the deviations from asymptotic 
form because it is exactly zero asymptotically: 



X D2(Mr7r)cOS 
1 
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(3{M^^yD2{M^ 



(12) 



Let us note, however that the expressions for intensity densities obtained above are based 
on leading twist expression for the amplitude (J^), in particular, the contribution of trans- 
versely polarised photon is neglected (it contributes at higher twist level). Since in the 
combination of intensity densities ([T2|) the leading twist contribution tends to cancel, the 
effect of higher twists (for which the cancellation is not expected) can be numerically 
large. 

The study of angular distribution allows to make separation^ of productions by lon- 
gitudinally (leading twist) and transversely (higher twist) polarised photons. Indeed, we 
can form combinations of intensity densities which get contributions only from 7^, so we 
minimise the contributions of higher twists. One of such combinations Q is, for example. 



d 



dM^ 



-(^0° + 75 1^2°) oc 
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X D2{M^^) cos 
1 



+ 



(13) 



The details of the derivation and the analysis of the data will be presented elsewhere. 



■^Under assumption of s— channel helicity conservation (SCHC), which holds with good accuracy ex- 
perimentally iQ. 

^Strictly speaking this combination gets also contribution from transverse polarisation of the photon, 
but it occurs only in i^-wave and therefore is expected to be small. 
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The advantage of the formahsm presented in this section is that it maximally uses 
the information on pion e.m. form factor F^(M^^) and pion phase shifts (5;(M^^) at 
low M^-TT which are very well known from low energy experiments. For example, from 
the known phase shifts 6i{Mt,t,) and Ss^Mt^^,) we conclude that the term proportional 
to cos[(5i(M7r7r) — Ss^Mj,.,,)] chaugcs the sign around M^-tt = 0.8 GeV, so to increase the 
sensitivity to this term it would be interesting to consider a "Mtttt asymmetry" for various 
observables: 

/ fO-8 l-Mma.^ \ d 

/ dM^^ - / c?M^^ -—— (an observable). (14) 

VJ2m^ Jo.8 ) aM^TT 

4. Expected values for the parameters 

By crossing relations the parameter is related to the third moment of the valence quark 
distribution in the pion 0: 

A2{(f) = iM^m = 1 f'dx {u^+ - u^+) . (15) 

D O JO 

If the parametrisation suggested in is used for the quark distributions in the pion, we 
obtain the values of the parameter A2{Q'^) listed in Table |I[ 



g2(GeV2) 


A2 


2 


0.110 


4 


0.099 


10 


0.089 


15 


0.085 



Table 1: Values of the parameter A2 as function of Q"^ , 
distribution in the pion suggested in ili 



obtained using eq. (15) and quark 



The value of Ai is constrained by the soft pion theorem [Q: 

A^{Q') = a^^\Q') - A2m , 

where a''2^ is the second Gegenbauer coefficient of the pion distribution amplitude. Unfor- 
tunately this parameter is not very well measured. In ref. [|13| the value of 02 = 0.19±0.13 
at Q2 = 2.4 GeV^ is quoted. 

Additionally the parameters y4i(Q^) and bi can be related to the second Gegenbauer 
coefficient a^2\Q'^) of the p meson distribution amplitude as 0: 

af{Q^) = A,{Q^)e'-^''^. (16) 

Up to now there is no direct experimental information about a2'\Q'^)- 
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The Mtttt dependence of the parameter 63 is estimated to be weak in the range of 
0.5 < Mtttt < 1-1 GeV and is then replaced by a constant 63. The parameters 6j were 
estimated in the instanton model of QCD vacuum |0 to be around the following values: 

= -0.75 GeV-2 
b2 = -0.75 GeV"^ (17) 
63 = 0.75 GeV^l 



5. Results of fits to the HERA data 

Recently data on di-pion mass distribution in hard exclusive reaction was measured at 
HERA by the reaction 

e + p ^ e + p + vr"*" + 7r~. (18) 

This section presents attempts to fit HERA data by the leading twist parametrisations 
described in the sections above. 

In order to study the evolution of the di-pion mass distribution three sets of data, 
the ZEUS photoproduction data (Q^ ^ 0) ig, the ZEUS low data (0.25 < < 0.85 
GeV^) [|T§ and the HI high data (2.5 < < go GeV^) |T§ are used. The mean W 
value for these samples is around 70 GeV, where W is the energy in the photon-proton cms 
(W'^ ~ Q'^/xBj — Q"^)- For the three samples, the di-pion mass distribution Q was studied 
in the ranges 0.55 < M^^ < 1.2 GeV, 0.55 < M^^ < 1.1 GeV and 0.50 < M^^ < 1.1 GeV, 
respectively. The main background to eq. (0) consists of events in which the proton 
diffractively dissociates into hadrons, and is estimated to be around 20 % for the ZEUS 
samples and around 10 % for the HI sample. This background does not distort the di-pion 
mass distribution discussed in previous sections if the mass of recoiled baryonic system 
is much smaller than the hard scale Q^. The contamination from the production of u 
(decaying into 7i~^7i~7i^) and (f) (decaying into K^K~ or 7r"'"7r~7r°) mesons was estimated 
to be few percents for the ZEUS samples and 7 % for the HI sample. This background is 
mainly situated at low mass M^^^^ < 0.6 GeV. 

Fig. presents the di-pion mass distribution (black points) for six Q"^ values: the first 
distribution corresponds to the ZEUS photoproduction sample, the two following ones to 
the ZEUS low sample and the three last ones to the HI sample. For the HI data, the 
distributions are corrected bin per bin for the production of uj and mesons. 

Firstly, we attempt to fit the HERA data with the simplest formula of the asymptotic 
form (^ where the di-pion mass distribution is directly proportional to the square of 
the pion e.m. form factor. This latter was recently precisely measured in low energy Q 
experiments (see also ||10|)) and we use directly these experimental data for \Ft,{Mt,t,)\. 
The result of the fits is superposed to the data in Fig. (black curves). The asymptotic 

^ Note that the di-pion mass distributions measured at HERA are not corrected for the transverse 
photon production. Nevertheless for greater than 2 GeV^ the longitudinal cross section dominates 
the transverse cross section. 

^ At leading twist the energy dependence factories out. The shape of the di-pion mass distribution is 
then energy independent in the present formalism. 
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form does not reproduce the distribution in photoproduction and at low Q'^ as expectedQ, 
but describes well the data at high Q^, the /ndf of the fits for the six bins of Fig. |l] 
being respectively 1625/25, 169/10, 69/10, 30.5/23, 16.6/21 and 11.6/18. We see that 
di-pion mass distributions in soft (low Q^) and hard (large Q^) regimes are essentially 
different and are described by different physics. Therefore the parametrisations designed 
for soft processes [0, |18|, |1^] are not relevant for hard processes, and the values of the 
parameters extracted from large data are not related directly to physical observables. 

Secondly, we study the deviations from the asymptotic form of the di-pion mass dis- 
tribution measured at HERA. As the formalism is valid in the hard regime ^ ^qcd' 
only the HI samples (Q^ > 2.5 GeV^) are considered. The parametrisation (p!0| ) is fitted 
to the data, with three free parameters: the normalisation N , and the parameters Ai and 
A2 . For the parameters &i , 62 , and 63 we use the values given by the calculation in the 
instanton model of QCD vacuum, see eq. (|1^. The result of the fits is presented in Fig. ^ 
as the black curves, and the values obtained for the parameters are listed in Table |^. The 
three fits have a /ndf value smaller than one. 



= 3.1 GeV2 


g2 = 7.2 GeV^ 


= 21.2 GeV^ 


= 1.59 ± 0.21 
Ai = 1.10 ± 0.18 
A2 = -0.34 ± 0.50 
X^/ndf = 20.6/21 


A^ = 1.11 ± 0.55 
Ai = 0.26 ± 0.49 
A2 = 0.22 ± 0.45 
xVndf = 16.5/19 


A^ = 0.66 ± 0.30 
Ai = -0.16 ± 0.34 
A2 = 0.18 ± 0.41 
X^/ndf = 11.3/16 



Table 2: Values of the normalisation N, the parameters A\ and A2, and the x^/ndf 
obtained from the fit of eq. ([T^ to the HI data [idiJ , using the values given in eq. ( \r^ ) for 
the parameters bi,b2, o-nd 63. 



We observe an indication for a decrease of the parameter Ai when increases, what 
is expected from QCD evolution, see eq. (^. In general the values of parameters Ai and 
A2 are in agreement with theoretical expectations (see section 4). However the errors of 
the parameters Ai and A2 are large, due to limited statistic available. Nevertheless the 
sensitivity of the data to the parameter Ai is encouraging. More precise data from the 
1997-99 years data taking at HERA would bring considerably more accurate information 
on these parameters which are related directly to the valence quark distribution in the 
pion, and pion and p meson distribution amplitudes. 

6. Discussion of the results and conclusions 

The analysis presented here is based on the leading order, twist-2 formalism. The power 
corrections (higher twist contributions) to the amplitude of the reaction (^ are systemat- 
ically neglected. The size of higher twist corrections might be rather large [^ , , so that 



''The asymptotic form (|^) is expected to work only in the hard regime ^ ^qcD' ^^^^ ^'^^ stress 
once more that the form of di-pion mass distribution of eq. (10) also makes sense only in the hard regime. 
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they constitute the main theoretical systematic uncertainty in our analysis. On general 
grounds the higher twist effects should be smaller for such observables as the shape of 
the di-pion mass distribution considered here, also the study of angular distribution of 
produced pions as discussed in section 3 can be used to minimise the contributions of 
higher twists. Another source of theoretical errors is that we fixed values of parameters 
bi by the model values (0). Hopefully more precise data will allow us to measure these 
parameters. 

To be on more safe side in respect to higher twist corrections we use the highest values 
of (Q^) = 21.2 GeV^ available in the data sample to evaluate observables related to 
partonic structure of the pion and p meson. The results for the third moment of valence 
quark distribution in the pion M3 (see eq. (1^)), the second Gegenbauer coefficients of the 
pion and p meson distribution amplitudes {a^^ and a'^^'' respectively) at {Q"^) = 21.2 GeV^ 
are presented in Table 



Quantity 


Our analysis 


Other sources 


Ms 


0.15 ± 0.35 


0.07 [12|, 0.085 ±0.005 


22 






0.02 ± 0.50 


0.14± 0.09 131 








-0.10± 0.20 









Table 3: Values of the third moment of valence quark distribution in the pion M3, and 
the second Gegenbauer coefficient of the pion and p meson distribution amplitudes (a^2\ 
af^ ) at (Q^) = 21.2 GeV"^ obtained in our analysis, and in 



IS, W 



Although, due to the limited statistic, the error bars for physical observables are large 
we see that the obtained values are in agreement with other experiments. This agreement 
is especially interesting because previously the restrictions for a\ and M3 were obtained 
from completely different measurements: the value of second Gegenbauer coefficient a\ 
was obtained in ref. [O from analysis of data on 771 transition form factor at large Q^, 



whereas the third moment M3 is restricted by vrA^ Drell-Yan data |T^, |2^0. We see that 
the formalism to extract the partonic structure of the pion and p mesons suggested here 
is complementary to already known methods. 

Our analysis allowed us to obtain for the first time an experimental estimate of the 
second Gegenbauer coefficient of the p meson distribution amplitude a'2' = —0.1 ± 0.2 
at = 21.2 GeV^. Unfortunately the precision of determination of ai'^ is still low 

to discriminate between different model predictions for this quantity (QCD sum rule: 
0.12 ± 0.06 m, 0.05 ± 0.01 [H and instanton model: -0.07 |, With increasing of 



statistics the accuracy of determination of the parameters related to the partonic structure 
of the mesons can be considerably improved. Also the studies of angular distributions 



^Lattice calculations of M3 are in agreement with results of analysis of 
0.10 ±0.02 H at ^ 21.2 



L2 231: 0.09 ± 0.03 24 and 
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of produced pions as discussed in section 3 can considerably increase the sensitivity to 
parameters of partonic structure of pions and their resonances. 

As a final remark we note that the formalism developed here can be easily gener- 
alised to the case of production of other mesons in hard exclusive reactions (e.g. KK, 
371, 4:71, KKnn, etc.). The study of meson mass and angular distributions of produced 
meson can provide us with rich information on partonic structure of these mesons and 
their resonances. Interesting prediction can be made 0, ^, ^ for the asymptotic shape of 
meson cluster mass distribution for the reaction: 

ll+P^ {hadrons) + N' , (19) 

where the mass of the hadron cluster M|- <^ Q^. In this case at asymptotically large 
we have: 



dN 
dMl 



oc R{Mx), with (20) 
a{e^e~ hadrons) 



It would be interesting to check experimentally this asymptotic formula and try to detect 
deviations from it. Generically, this deviation is described by the vacuum correlator of 
two light-ray opertors. We showed here that the asymtotic expression (^) works pretty 
well at low Mx (for > 1 GeV^) where the 27r channel dominates over all other hadronic 
channels, it would be interesting to check how good this formula works at higher Mx- Also 
it would be interesting to see whether the formula like (^) is applied also for multimeson 
partial cross section. 
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Figure 1: Di-pion mass distribution for six values. The black points correspond to 
measurements from the ZEUS [7^ and HI j\r^ experiments. The curves show the 
result of the fit of the asymptotic form 
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Figure 2: Di-pion mass distribution for three values. The black points correspond to 
measurements from the HI /jT^ experiment and the curves show the result of the fit of 
eq. 
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